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Abstract 

We consider joint Tikhonov- and Lavrentiev-regularization of control 
problems with pointwise control- and state-constraints. We derive er- 
ror estimates for the error which is introduced by the Tikhonov regular- 
ization. With the help of this results we show, that if the solution of 
the unconstrained problem has no active constraints, the same holds for 
the Tikhonov-regularized solution if the regularization parameter is small 
enough and a certain source condition is fulfilled. 
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1 Introduction 

In this paper we consider problems that can be interpreted either as optimal 
control problems or identification problems (inverse problems). Let D' C -D be 
bounded domains in (iV = 2, 3), J7 = L'^{D) and F be a Hilbert space. For 
a compact and linear mapping S : U Y , an element yd & Y and bounded 
measurable functions 6 : _D — > R and : Z?' — > R we consider the constrained 
minimization problem 

min||5u — ydip s.t. < u < 6 a.e. on £>, S'm < t/j. a.e. on D' (P) 

The space Y is called data space in inverse problems and state space for control 
problems. Moreover U is considered as solutions space or control space. Our 
specific feature is the presence of two types of different inequality constraints. 
The first one 

Q < u <h a.e. on D 
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describes a set of physically meaningful solutions in inverse problems or a set of 
admissible controls. The second one 



Su < ip. a.e. on D' 

ensures that the data (the state) is pointwise in a reasonable range. In the 
inverse-problem context, the constraints on u and Su model prior knowledge on 
the solution and the measured data, respectively and hence, they shall lead to 
more reliable reconstructions. 

Before we start to analyze this problem, we will give a specific example. Let 
S : u i—f y he the solution operator of the boundary value problem 

—Ay = u on D, y = on dD 

with a Lipschitz domain D. We specify the spaces U ^ Y = L'^{D). The set 
D' C D can be for instance an inner subdomain. However, it would also be 
possible to discuss sets D' containing only finitely many points. Of course our 
setting fits also to more challenging problems. 

Next, we will discuss the regularization of such a problem. Let us first men- 
tion that the introduced problem (|P| can be well-posed or ill-posed. One can 
easily construct situations, where the set of functions u satisfying the inequality 
constraints is empty or consists of exactly one point. However, such situations 
are not in our focus. We will assume later that the set of feasible functions 
u G U has an inner point with respect to the L°°-topology. 

Since the set of feasible functions u is weakly compact in [/, the existence of 
at least one function u satisfying the inequality constraints ensures the existence 
of a solution of ([£1 . Consequently, existence of solutions of will not be an 
important issue in this paper. 

Let us now think about uniqueness of solutions. Using standard arguments 
one can show the uniqueness of the state (resp. data) y := Su. However, the 
uniqueness of u is only guaranteed if S is injective. Of course, a standard 
Tikhonov regularization guarantees uniqueness of solutions of the regularized 
problems. Therefore, the discussion of uniqueness aspects seems to be com- 
plete. However, there is another uniqueness aspect occurring even in the case of 
finite dimensional spaces. Solution of minimization problems are analyzed and 
computed by means of optimality conditions. The most convenient form of op- 
timality conditions includes Lagrange multipliers, i.e., dual or adjoint variables. 
It may happen that there exist a subset D" C D' where more than one inequal- 
ity holds as equality. Typically this effect is connected with nonuniqueness of 
the dual variables on the corresponding set D" . For the construction of specific 
examples we refer to [5]. Let us mention that Lagrange multipliers associated 
with the inequality Su < ip are in general only Borel measures, see Casas [I]. 
This low regularity of the dual variables was the motivation in [lOj to introduce 
a second Lavrentiev type regularization. 

Motivated by this argumentation we will study in this paper a family of 
Tikhonov and Lavrentiev regularized problems: 

min \\Su — yd\\^ + s.t. < u < b, a.e. on D, Xu + Su < ip a.e. on D' 

with a, A > 0. We denote with u a solution of and with u^ a solution 



of ( I . Note that the solution of ( ) is always unique if it exists 



2 



Next, one can study a lot of different regularization errors. People coming 
from inverse problems usually study the solution behavior for a | 0. Since the 
inverse theory heavily uses representation formulas, pointwise inequalities are 
undesired. 

People coming from optimal control state that the Tikhonov term represents 
the costs for the control and assume that the parameter a is a given, fixed model 
parameter. Hence, it would be enough to study the behavior for A | for a 
fixed Tikhonov parameter a > 0. However, even in the optimal-control context 
one may take the position that the problem (jP]) is the one which shall be solved 
and that the Tikhonov regularization is employed only to stabilize the problem. 

Of course, error estimates for the whole regularization process are of high 
interest. However, this problem seems to be challenging. In this paper, we will 
only contribute a little bit in answering this question. 

Let us summarize our plans. The problem (jP| is the one which we want to 
solve. Therefore, we are interested in estimates for the error m — m^. We split 
this error into the Tikhonov error and the Lavrentiev error: 

\\u~u^j< \\u-vu + • 

Tikhonov error Lavrcnticv error 

Estimates for the Lavrentiev error can be found in [Tni[51|21[n]- Here, we focus on 
estimates of the Tikhonov error. The analysis of this error is well developed in 
the framework on inverse problems for problems without inequality constraints. 
However, for the constrained case there are only few results ^5, Section 5.4]. 

The paper is structured as follows: In Section 2 we will state some prelim- 
inary results. Error estimates for the constrained Tikhonov regularization are 
located in Section 3. The activity of inequality constraints is analyzed in Section 
4. Error estimates for the Lavrentiev regularization can be found in Section 5. 
The verification of the assumptions for a distributed elliptic control problem 
and a Fredholm integral operator is presented in Section 6. 

2 Preliminary results 

First, we introduce notations for the admissible sets for and 

U^d^{ueU \0<u<b, Su<il)}, 

^{ueU\0<u<b, Xu + Su<ij}. 

Note that, due to < u we have J/ad'*' C C/ad^ for A > /i. With this notation we 
reformulate 

min - s.t. w G C/ad (P) 

mm\\Su~yaf + a\\u\\^ s.t. ueU^d^. (P^) 

Now we state some preliminary results: 

Lemma 2.1. Let the operator S be linear and continuous. Then the sets J7ad 
and C/ad^ are closed, convex and bounded. 



Pi 
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Definition 2.2. The constrained pseudo inverse of an operator S with respect 
to a convex and closed set C is defined via 

\\Sciyd)\\ =min{||M|| | u e a.rgmm\\ Su - yd\f}. 

c 

In other words, S'^{yd) is the minimizing element of the residuum which has 
minimal norm. 

The operator is a nonlinear operator with the following properties. 
Lemma 2.3. If C is non-empty, convex, closed and bounded, then 

DiS+) = Y. 

Proof. Note that C is weakly sequentially closed and hence a minimizing se- 
quence of \\Su — j/dlP in U has a weak accumulation point. By lower semiconti- 
nuity we see that this accumulation point is indeed a minimizer, and hence, the 
constrained pseudo inverse exists. □ 

The following well known proposition shows continuity in the Tikhonov pa- 
rameter a: 

Proposition 2.4. depends continuously on a for a > 0. 

Proof. We drop the superscript A because it is fixe d he re. Consider Ua and Ufs 
for a,/3 > 0. Since A is fixed, up is admissible for (P^) and vice versa. Hence, 
we can insert these elements in the corresponding variational inequalities and 
obtain 

{S* Sua + aua — S*yd, up — u^) > 
{S*Sup + f3up - S*yd, Ua - Up) > 0. 

We get 

{S*SUa + CeUa — S*Sup — I3up, Up — Ua) > 

and obtain 

\\SUa - SupW'^ < {aUa - f3up,Up - Ua) < \a ~ l3\\\Ua\\\\up - Ua\\ ^ l3\\up - Maj^. 

Since the left hand side is positive we get 

\\up-Ua\\ < ^ ll^aH 

and the right hand side converges to for /? ^ a. □ 



Finally, we state a characterization of the solution of ( P^ ] by means of a 
projection formula: 



Lemma 2.5 (Projection formula). The solution u^ of (Pa) is characterized by 
Ua^Pu.,^4 ■ (1) 
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Proof. We rewrite {P^ ) with the help of the indicator function Ijj ,>■ as 



minll^w - yd\\^ + a\\u\\^ + Ijj^^x{u). 



With the help of subgradient calculus we get from optimality of 

Since d{\\ ■ |p + — Pu^^>- (see, e.g. [H]) this gives 

S*{Sui-yd) 



□ 



3 Error estimates for constrained Tikhonov reg- 
ular ization 

Now we establish error estimates for the term — Resembling results can 
be found in [5] . Due to the structural differences between inverse problems and 
optimal control we state the result in our terminology with an explicit constant 
and present a full proof. 

Theorem 3.1. Let u he a solution of [P\) and denote with Pjj^d projection 
onto J7ad- Moreover, let the following source condition be fulfilled: 

3w€Y : w = Pu^^iS*w). 

Then it holds: 

\K-u\\ < ^\\w\\ + ^\\Su~y,\\ (2) 
\\Sul-yd\\ < 2a\\w\\ + \\Su^yd\\ (3) 

Proof. By definition we have € C/ad and u € C/ad- Hence, by optimality of 
we have 

\\Sul-y,r+a\Kr<\\Su-y,r + a\\ur. 
Rearranging yields 

\\Su"^~ydf + a{\\u'if-\\uf)<\\Su-y,f 

which we extend to 

|j5<-y,|r+a(||<f-|1^2|p-2(5*«;,<-^2)+2(5*«;,<-u))<||5u-y,f. 

(4) 

Since u = Pu^dS*w we have for all u e C/ad 

{S*w — u, u — u) < 0. 
Using this with m = we get 
IKIP - Wuf 2(S*w, < - > IKIP - - 2{u, < - u) = IK - uf. 
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We further estimate from ^ 

\\Sul - ydf + -u\\^ + 2{S*w, < - u)) <\\Su- yaW^ 

^ ||5<-2/d||' + 2(au;,5<-y<i)+a||<-u||2 < \\Su - ydf + 2{aw, Su - yd). 
Completing the squares by adding ||q!1(;|P gives 

WSu'i-yd + awf +a\\u'i-uf < \\Su - ya + awW^ . 
On the one hand this leads to 

\K-u\\'<i\\Su~yd + aw\\' 
which gives by taking square roots and using the triangle inequality 

\K-u\\<^\\Su-yd\\+V^\\w\\. 
On the other hand we conclude 

\\Sul-yd + awf < \\Su-yd + aw\\^ 

which implies 

\\Sul-yd\\ < \\Su-yd\\+2a\\w\\. 

□ 

The estimate ^ motivates the following parameter choice: If yd is not in 
the range, we see that the right hand side in ^ is smallest for 



and hence, is a reasonable choice for the choice of the Tikhonov parameter a if 
the quantities were known. Nonetheless, the error estimate ^ is useful for the 
determination of the total error — 



4 Activity of the constraints 

In this section we investigate the following situation: Assume that the optimal 
solution u of jPj has no active inequality constraints. That means that in fact 
we would get the same solution without imposing any inequality constraints. 
However, the formulation with additional inequality constraints is reasonable 
since a solution of an unconstrained inverse problem can violate these constraints 
for noisy data. 

Now the question arises: Can we expect a solution without active constraints 
also for the purely Tikhonov-regularized problem for small regularization param- 
eters? In fact this can be shown with the help of the estimates of Theorem 13. II 

Theorem 4.1. Let S : U ^ L'^{D) ^ Y with dense range and assume that 

Un u in U Sun — > Su in L°"{D') (5) 
yn^y mY ^ S*y„^Sy inL°°{D) (6) 

Let u be a solution of lP\) such that for some t > it holds that t < u < b — t 
and Su < ip ~ t hold. Moreover, let there be w ^Y such that u ~ S*w. 

Then there exists > such that for every a < the solution of (P^) 
also fulfills T < < b ~ T and Su^ < ijj ~ t . 
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Proof. Since u does not have active constraints and S has dense range, it holds 
Su = Ud- To see this, note that u fulfills 

{S* {Su -yd),u-u) >0 for all u e f/ad 

and let us assume that there exists v G U but v ^ J7ad, such that 

{S*{Su~yd),v-u) <0. (7) 

We may approximate v by bounded functions and hence, there exists a bounded 

V such that 

{S*iSu-yd),v-u) <0. 

Now observe that for 6* > small enough, the function — 9v + {1 — 9)u is in 
C/ad since u does not have active constraints. We see 

{S*{Su-yd),v' -u)^0{S*iSu-yd),v-u) <0 

which contradicts the optimality of u for Hence, ([7]) has to be fulfilled for 
all w e [/ and this shows that u is also a solution of the unconstrained problem. 
Since S has dense range, this optimal value for this problem is which shows 
Su = yd- 

Moreover, we have by assumption u = S*w = Pu^^S*w. We conclude from 
Theorem [XT] that 

ll<-^ll<x/^IIHI- 

This implies u'^ —> u in and by assumption we know that 5*^^ — > Su in 
L°°{D'). Because of Su < ip — t, we find Su'^ < tjj — t/2 for sufficiently small 
a. Hence, we can use the projection formula ([T]) with P[o,b] instead of Pjj^d fo^' 
small a and get 

<=P[o,6](-5*i(5<-yrf)). 
Again from Theorem 13.11 we find 

\\Sul-yd\\<2a\\w\\. 

Consequently, ||(5'u" — yd)/ce\\ is uniformly bounded in Y. Let us now take 
an arbitrary sequence {a„} with q;„ ^ for n to oo. Since {Su'^^ — yd)/ctn 
is uniformly bounded in we can find a weakly convergent subsequence in 

Y denoted by the index n' . The assumption on S* yields strong convergence 
of u° , in L°°{D). We already know that , converges in U — L'^{D) with 
limit u. By a standard argumentation we obtain that u*^ —^ u in L°°[D) for 
a — > 0. Consequently, the control constraints are also inactive for sufficiently 
small a. □ 

Let us note that the first argumentation in the proof works also for the 
weaker assumption that m„ — > u in LP{D) implies Sun — s- Su in L°°{D') for a 
sufficiently large p. However, the assumptions for the adjoint operator cannot 
be weakened. Thus, the practical benefit of that generalization is only small. 

The above result is also of interest in the case of ill-posed problems. Here 
we may ask the question if additional constraints on the solution in the mini- 
mization of the Tikhonov functional may destroy the optimal convergence rate. 
As we will see below, this is not the case in our setting. Similar to Theorem 13. II 
we can state the following result: 



7 



Theorem 4.2. Let u be a solution of IP]) with Su — yd and let he such that 
\\v^ ~ UdW < 5. Moreover, let the following source condition be fulfilled: 



Then it holds for 



that 



3w^Y : Pu,^iS*w). 



■"^ - argmin||S'ii- + Q!||w|l2 (8) 



\\Sui-yd\\ < 2a\\w\\+5 

It also holds that all constraints get inactive for small a if the true solution 
of noisy data does not have active constraints: 

Theorem 4.3. Let S : U ^ Y fulfill the conditions and ([^. Let u be 

a solution of lP\) such that for some t > it holds that t < u < b — t and 
Su < ip — T holds. Moreover, let there he w (zY such that u = S*w and let 
be defined by Finally let a{S) be a parameter choice rule such that 

a(S) 0, — -:- ^0 for 6^0. 
a[d) 

Then there exists Sq > such that for 5 < So it holds that r < u^^^^ < b — t and 

^'^a{5) < "4^ ^ T holds. 

Proof. Due to the parameter choice we get w^j-i-) u strongly in for S ^ 0. 
Similar to the proof of Theorem 14.11 we conclude that the state constraints are 
not active for sufficiently small S and hence, the projection formula 

holds. Now the claim follows similarly to Theorem 14. II □ 

Finally we state the following converse result which shows that the source 
condition u = Pu^^S*w follows from weak convergence of the regularized solu- 
tions together with a mild decay of the discrepancy. 

Theorem 4.4. Let u'^ ^ u and \\Su^ — yd\\ — 0{a) for a ^ 0. Moreover, 
we assume that S fulfills |^ and Then there is a function w such that 

Pu^iS*w = u. 

Proof. Since ||S'u° ~ yd\\ = 0{a) we have that 

\\Sul-yd\\ , 
hmsup < oo. 

Hence there is a sequence a„ — > and an element w such that 

Sua^ - yd 
^ —w. 

Oin 
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By the projection formula ([T]) we have 




S*{Su\ 







yd) 



) 



= u, 



,0 



'a. 



For the right hand side converges weakly by assumption. With a discussion 
similar to that one in Theorem 4.1 we obtain the strong convergence of the left 
hand side of the last equation. Since the both sides converge weakly we have 
by uniqueness of the weak limit 



5 Lavrentiev regularization 

We discuss now the additional Lavrentiev regularization. The motivation of this 
second regularization is to improve the properties of the adjoint problem. In 
this section we will investigate two different situations. In the first part we will 
sketch the estimation of the Lavrentiev error in the general case. In the second 
one we will again discuss the situation of Theorem 14. II Then we will be able to 
show better convergence results. 

Let us start with the general case. The discussion for fixed a > and a 
specific problem can be found in [3]. Next, we reformulate the assumptions of 
that paper for our more general setting. 



The assumption (llOp means that there exists a Slater point with respect to 
the state constraints. In the second part of this section we will benefit from 
the stronger assumption that u itself has this Slater property. In Section 2 we 
already mentioned that the set t/ad^ of admissible u becomes smaller when A 
becomes larger. However, one can show the existence of at least one admissible 

control for A < n-n . 

The error estimates are obtained by the following technique: 

1. Write down the necessary optimality conditions for vP^ and as varia- 
tional inequalities. 

2. Take as test function in the optimality condition of vP^. 

3. Construct a convex linear combination = au + (1 — cr)u^ which belongs 
to C/ad^- Take this function as test function in the optimality condition of 
<■ 

4. Add both inequalities and estimate all terms. 

The resulting error estimate can be found in [31 Theorem 5.4]: 



Pu^aS*w = u. 



□ 



S -.U = L^{D) ^ L°^{D') continuously. 
There exists ii E U,t > such that < u < b, Su < tjj — t. 



(9) 
(10) 
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Lemma 5.1. Let the assumptions (0) and I110\) be fulfilled. Then there ex- 
ists a constant c > such that for X < „.,, the error of the Lavrentiev 

regularization can be estimated by 

\\ul-ui\\<c-. 

a 

This general result has an essential drawback: If a becomes small, then A 
has to be very small to ensure a certain accuracy. 

Let us now assume that u itself has the Slater property 

Su < ip — T. 

In contrast to Section 4 we do not require an inner point property with respect 
to the control constraints. 

Theorem 5.2. Let S fulfill assumption ^ and let u be a solution of fP)) such 
that the Slater condition Su < ip — t and a source condition u — Pu.^^ {S*w) are 
fulfilled. Then it holds for sufficiently small A and a that 



Proof. From Theorem 13. II we know the convergence of vP^to umU . Using the 
mapping property S : U = L^{D) — > L°°{D'), we can show as in the proof of 
of Theorem rO that 

Sul<ij- t/2. 

holds for sufficiently small a. From Lemma 15.11 we know that tends to 
for A ^ 0. Due to the properties of S we find 

Su^<^ 

for sufficiently small A. For A < ^ we get 

+ Xul < ^. 

Consequently, m° is feasible for the problem {P^) and is feasible for the 
problem (P^) for a and A small enough. Testing the optimality condition for 
{Pa) and with uj^ and respectively, yields, similar to the proof of 

Proposition EH that = u° . □ 

Let us remark that the Lavrentiev regularization is also used with different 
sign, i.e., 

Su — Xu < ip. 

Then, the set of admissible controls becomes larger for larger A. Hence, the 
Slater condition ensures the existence of feasible controls for arbitrary positive 
A. In this case we have no smallness condition for A. The general result for the 
Lavrentiev regularization can be found in \T, Theorem 3.3]. The discussion for 
the specific case of Theorem 15.21 can be done completely analogue. 

Both approaches have their specific advantages. As we have already seen, the 
first approach (plus sign) yields solutions that are feasible for the unregularized 
minimization problem. However, we have to deal with a smallness condition for 
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the Lavrentiev parameter A. The second approach (minus sign) does not need 
an additional smaUness condition, but the solutions are in general not feasible. 

Finally we combine the results of Theorem 13.11 Theorem 14.11 and Theo- 
rem 15.21 to obtain an estimate for the total error of joint Tikhonov-Lavrentiev 
regularization under strong assumptions: 

Theorem 5.3. Let S fulfills the assumptions 0) and (0j and let u he a solution 
of IP]] such that there exists r > such that t < u < b ~ t, Su < ~ t and a 
w £ Y such that u — S*w. Then it holds for A small enough that the solution 



< of iP^) fulfills 



=^0{y/a) for a^O. 



Proof. We split the total error as 

\\u~u^J < \\u~ul\\ + \\ul-u^^\ 



and observe that due to Theorem 15.21 the second term vanishes if A is small 
enough and that the first term behaves like 0{y/a) due to Theorem 13.11 and 
Theorem O □ 



6 Verification of the assumption in special cases 
6.1 An elliptic control problem 

In this section we will discuss the example from the introduction: Let S : u y 
be the solution operator of the boundary value problem 

—Ay = u on D, y = on dD 

with a Lipschitz domain D C R^, N E {2, 3}. We specify the spaces U = Y = 
L^{D). The set D' C D is assumed to be an inner subdomain. Let us now 
check the assumptions: 

1. By the Lax-Milgram Lemma we obtain easily the existence of a unique 
solution y e Hq{D) for every u £ U. Due to the embedding Hq{D) ^ Y, 
the operator S is well defined 

2. Let us first mention that the operator S is selfadjoint for our specific choice 
of spaces. For the uniform boundedness of solutions of the elliptic problem 
we refer to Stampacchia [T^ . 

3. The mapping property in Theorem 14 . 1 1 that y„ y in K implies S'*y„ — > 
S*y in L°°{D) can be obtained by the following argumentation. Weak 
convergence yn ^ y vaY = L^{D) implies strong convergence in yn —> y 
in W-^^P{D) for 2 < p < oo for iV = 2 and 2 < p < 6 for iV = 3. Now the 
desired result follows again from [12, Theorem 4.2]. 

4. The operator S has dense range in Y. That can be verified by the following 
argumentation. The space C^{D) is dense in Y. Moreover, all these 
functions belong to the image of S. 

5. The Slater condition can only be checked (analytically or numerically) if 
one specifies the data. The Slater property for u is like the source condition 
an a priori assumption for the solution. 
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Consequently, we can apply all results of our paper to that example. Only 
the Slater property and the source condition are a priori assumptions. All other 
assumptions of that paper are satisfied for our example. 

6.2 A Fredholm integral operator 

Another class of examples in which the properties 1.-4. are fulfilled and which 
models several inverse problems is given by Fredholm integral operators 7J or 
[1 Chapter VI.]. Let £) be a bounded domain in and \etU = Y = L^{D). 
For a Lipschitz continuous function k : D x D ^ H we consider Su{x) — 
k{x, x')u{x')dx' . Let us check the assumptions: 

1. Since D is bounded, continuity of k implies that S : L'^{D) — > L^{D) 
compactly. 

2. The uniform convergence of 5w„ to Su follows from 7i„ ^ u by 

||S'u„ — 5w||oo < sup / \k{x^x')\\un{x') — u{x')\dx' 

x£D J D 

< sup \k{x,y)\y/\D\\\un- u\\2. 

3. Uniform convergence of S*yn to S*y follows from y„ y by observing 
that \\yn — 2/II2 is bounded and that the mapping a: 1-^ k(x' , x)y{x')dx' 
is Lipschitz continuous for every y G LF'{D). 

4. Since 5* is compact is has a representation via its singular value decom- 
position Su = '^.fi<yn{uTipn)4'n with uou-uegative singular values which 
decay to zero. We see that the range of S is dense if S does not have zero 
as a singular value (i.e. it is injective). 

Again, the Slater condition and the source condition are a priori assumption 
for the solution. 

Conclusion 

In this paper we studied the simultaneous Tikhonov and Lavrentiev regular- 
ization of an optimal control problem with control and state constraints. We 
derived error estimates in the general case and the main tool was a source con- 
dition which resembles the classical one used in the inverse-problem context. 
With the help of this error estimate we could prove, under additional assump- 
tions, that for sufficiently small Tikhonov parameter the Tikhonov- regularized 
solution does not have active inequality constraint if the original solution has the 
same property. Moreover, it was shown that for the Tikhonov- and Lavrentiev- 
regularized problem the solution coincides with the Tikhonov-regularized solu- 
tion for sufficiently small regularization parameters if the unregularized solution 
has a Slater property. One may conclude that using additional physically moti- 
vated inequality constraints in the context of the regularization of inverse prob- 
lems is a good idea for two reasons: For larger regularization parameters they 
may yield reconstructions that are more meaningful and asymptotically they do 
not destroy optimal convergence rates (and even become inactive if they will 
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be in the limit). On the other hand, additional inequality constraints lead to 
minimization problem which may be harder to solve. However, there are pow- 
erful algorithms available which allow a numerical treatment of such problems 

(cf. [HE]). 

Moreover, our results provided a little insight in the problem of error esti- 
mates for joint Tikhonov-Lavrentiev regularization for optimal control problems 
and obtained preliminary estimates. Further research on the interpretation of 
the source condition and the Slater condition in particular contexts seems nec- 
essary. Finally we showed that the main assumptions for our setting are fulfilled 
for a distributed elliptic control problem and for some Fredholm equations of 
the first kind. 
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